arXiv: 1501.04260v 1 [cs.SY] 18 Jan 2015 
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Abstract —In this paper we study disease spread over a 
randomly switched network, which is modeled by a stochas¬ 
tic switched differential equation based on the so called 
A'-intertwined model for disease spread over static networks. 
Assuming that all the edges of the network are independently 
switched, we present sufficient conditions for the convergence 
of infection probability to zero. Though the stability theory 
for switched linear systems can naively derive a necessary and 
sufficient condition for the convergence, the condition cannot 
be used for large-scale networks because, for a network with 
n agents, it requires computing the maximum real eigenvalue 
of a matrix of size exponential in n. On the other hand, our 
conditions that are based also on the spectral theory of random 
matrices can be checked by computing the maximum real 
eigenvalue of a matrix of size exactly n. 

I. Introduction 

Network epidemiology, a branch of mathematical epi¬ 
demiology which aims to mathematically understand dis¬ 
ease spread over networks, has been attracting ever-growing 
attention [1], One of the reasons comes from the finding 
that the heterogeneous structure of real world networks, 
such as an uneven distribution of the number of neighbors 
of individuals, cannot be ignored to perform an accurate 
prediction of epidemics [2], Another reason comes from 
recent advancements of the technologies for collecting a 
massive amount of data from real human interactions, which 
enable researchers to build accurate mathematical models of 
human networks (see [3] and references therein). One of the 
recent achievements of network epidemiology, for example, 
is the effective prediction of the pandemic of H1N1 influenza 
on 2009 [4], 

In the last decade, the major emphasis of network epidemi¬ 
ology has been on the disease spread over static networks [1], 
This time-invariance assumption enables us to model disease 
spread as a time-homogeneous Markov process, for which an 
effective approximation by a constant-coefficient differential 
equation, called an A'-intertwined model, is available [2], 
[5]. For example it is known [2], [5] that the evolution of a 
certain epidemic model over a static network is characterized 
by only the maximum real part of the eigenvalues of the 
adjacency matrix of the network. 

However, recently, several simulation studies have re¬ 
ported that the dynamical nature of human interactions can 
greatly affect the way epidemics evolve and that static 
approximations of such dynamical networks may result in 
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inaccurate conclusions [6]. For example, Vazquez et al. [7] 
observe that the speed of disease spread can be substantially 
slowed down on randomly switched networks than on static 
ones. The result by Volz and Meyers [8] suggests that the 
frequency of the change of the topology of network structure 
is an important factor for the emergence of endemic. In 
[9], it is numerically confirmed that so-called memory effect 
of networks can slow down disease spread. Despite those 
simulation results, rigorous understanding of disease spread 
over dynamical networks based on mathematical analysis is 
very limited. 

This paper gives an analysis of disease spread over 
randomly switched networks with the stability theory for 
switched linear systems called Markov jump linear sys¬ 
tems [10] and the spectral theory of random graphs (see, 
e.g., [11]). We provide sufficient conditions for epidemics 
to eventually die out over randomly switched networks 
having independent and Markovian edges. This paper can 
be considered as a stochastic counter-part of the recently 
submitted paper [12], where the authors study epidemics over 
deterministically switched networks. 

One of the major difficulties of the problem is in the 
computational cost. The stability theory for Markov jump 
linear systems in fact gives a necessary and sufficient con¬ 
dition for epidemics eventually dying out as we will see 
later (Proposition 3.1). Unfortunately, checking the condition 
requires finding the maximum real eigenvalue of a matrix 
with size 2"(" _I '>/ 2 , where n denotes the number of the agents 
in a network. On the other hand, the size of the matrix 
appearing in the proposed conditions is equal to n, which 
allows us to apply the conditions for large-scale networks. 

This paper is organized as follows. After giving the 
notation used in this paper, in Section II we give the model 
of disease spread over randomly switched networks. Then 
Section III gives sufficient conditions for epidemics dying out 
over randomly switched networks having independent edges. 
Examples are presented in Section IV. Finally in Section V 
we extend the obtained results to the case when networks 
are modeled by weighted graphs. 

A. Mathematical Preliminaries 

We let /„ denote the n x n identity matrix. The subscript n 
will be omitted when it is obvious from the context. A real 
vector v is said to be nonnegative, written v > 0, if it has only 
nonnegative entries. By t n and l„ m we denote the u-vector 
and the n x m matrix whose entries are all one. A square 
matrix is said to be Metzler if its off-diagonal entries are 
nonnegative. We denote the Kronecker product of matrices A 
and B by A <g> B. 


The Euclidean norm of x £ R" is denoted by ||x||. Also we 
define the 1-norm of xby ||x||i = £" =l |x, |. The maximum real 
part of the eigenvalues of A is denoted by 77 (A). We say that 
A is Hurwitz stable if 77 (A) < 0. The matrix measure [13] 
of A £ R nx " is defined by /i(A) = lim/^o (||/+A/z|| — 1 )/h, 
where ||-|| denotes the maximum singular value. When A is 
symmetric and thus has only real eigenvalues, its maximum 
eigenvalue is denoted by A max (A). 

An undirected graph is a pair Sf = ('V. S’), where "V is 
a finite and nonempty set and S is a set consisting of 
distinct and unordered pairs (i,j) of i. j £ “V. We call the 
elements of "V vertexes and those of & edges. The adjacency 
matrix A^ of Sf is defined as the square {0, l}-matrix of 
size \Y\ whose (i,j) entry is one if and only if (i,j) £ S. 
Clearly the adjacency matrix of an undirected graph is 
symmetric and has zero diagonals. We say that j £ 'V is 
a neighbor of i £ f' (or ; and j are adjacent) if (i,j) £ V. 
The degree of i is defined as the number of the neighbors 
of i. 

The expectation of a random variable is denoted by £[•]. 
All the Markov processes appearing in this paper are as¬ 
sumed to have finite state spaces. 

II. Epidemic Model over Randomly Switched 
Networks 

In this section we give a model for the spread of disease 
over randomly switched networks using switched differen¬ 
tial equations. The model is based on a model called TV- 
intertwined model originally proposed for epidemics over 
static networks [5]. Then we will review some results on the 
stability of switched linear systems that will be used in this 
paper. 

A. TV -intertwined Model for Static Networks 

We first give an overview of the TV-intertwined model [5] 
for so-called susceptible-infected-susceptible types of dis¬ 
ease. Let a population be modeled by an undirected graph Sf 
over vertexes {l,...,n}. We regard each vertex as an agent 
who can be infected and also can transmit disease to its 
neighbors. In the model, at each time t > 0, each vertex 
can be in one of the two states: susceptible or infected. We 
assume that, when a vertex i is infected, the transition to the 
susceptible state occurs following a Poisson process with 
rate 5,-, called curing rate. On the other hand, if i is in the 
susceptible state and one if its neighbors is in the infected 
state, then i will make a transition to the infected state by 
following a Poisson process with rate If, called infection 
rate. All the Poisson processes are assumed to be indepen¬ 
dent. Throughout this paper we assume that population has 
homogeneous infection rate and curing rate S, i.e., we 
assume that /3i = • • • = j3„ = j3 and 5i = • • • = 8 n = 8. 

Let pft) denote the probability that the vertex i is infected. 
Define p = [p\ ■■■ p„\ T and P = diag(/?i,... ,p n ). Then the 
TV-intertwined model [5] reads 

Sl)p- fiPAyp, ( 1 ) 


where p( 0) = pq£ [0,1]". We notice that clearly p(t) > 0 for 
every t > 0 . 

B. Switched TV -intertwined Model 

Based on the TV-intertwined model (1) we can readily 
state our model of disease spread over randomly switched 
networks. We assume that our randomly switched network is 
modeled as := {^a(t)}t>t), where a = {d(f)};>o is a time- 
homogeneous Markov process taking its values in {1,... ,7V} 
and ,..., Sfjv are undirected graphs over vertexes {1 ,... ,n}. 
Examples of such random switched graphs include the activ¬ 
ity driven networks [14], graphs with edge swapping [ 8 ], and 
temporal exponential random graphs [15]. Then we model 
the disease spread over the randomly switched network 
by the switched TV -intertwined model 

E: % = ^ A 'U) ~ 5I ^ P ~ P pA %nP’ 

where p{ 0) = po £ [0,1]" and cr(0) = Oq £ {1,...,7V} are 
arbitrary constants. 

The principal aim of this paper is to give conditions under 
which the zero equilibrium of E is stable, i.e., the infection 
probability p converges or stays close to the origin. We 
introduce the following definitions. 

Definition 2.1: We say that E is 

1) mean stable if there exist C > 0 and e > 0 such that, for 
every po and Co, it holds that £’[||/?(f)||] < Ce“ w ||po||; 

2) almost surely stable if, for every po and do, it holds 
that P(lim,^ OT ||p(f)|| = 0) = 1. 

The following linearized model 

= (/U ^> ~ 5I)P (2) 

of E plays as an important role as the one for static 
networks does in [16], [17], where the authors propose 
the optimal vaccination strategy for the disease spread over 
static networks. We let the initial conditions of E given 
by constants 7?(0) = po £ [0,1]" and <r(0) = do £ {1,... ,7V}. 
The next lemma is easy to see but fundamental. 

Lemma 2.2: If po = po and do = do, then we have 
11/7(011 1 < ll/KOIIi for every / > 0 with probability one. 

Proof: Let e(t) = p(t) — p(t). Then e(0) =0 and also 
de/dt = /3 PAcg^p > 0 because all of /3, P, Ag^, and p are 
nonnegative. Therefore p(t) < p(t) entry-wise for every t > 
0. Finally multiply ij to the obtained inequality from the 
left to obtain the desired inequality. ■ 

Therefore, the study of the stability of E can be reduced 
to that of the switched linear system E, for which various 
results from the theory of switched linear systems [ 10 ] are 
available. We will review some of the results in the next 
section. 

C. Stability of Markov Jump Linear Systems 

As in the previous section let d = {d(f)} r >o be a time- 
homogeneous Markov process taking its values in {1,... ,TV} 


and let Ai,...,Aw £ R" x ”. Consider the switched linear 
system (called Markov jump linear system [10]) 

dx 

dt — A cs(t)X, (3) 

where x( 0 ) =xo £ K" and cr( 0 ) = Oq are arbitrary constants. 
The stability of Markov jump linear systems is defined in 
the following standard way. 

Definition 2.3: We say that (3) is 

1) mean stable if there exist C > 0 and e > 0 such that 
£[||x(f)||] < Ce~ e, ||xo|| for all xo and Ooi 

2) almost surely stable if, for all xo and Ob, there holds 
P(lim r ^oo||x(f)|| =0) = 1. 

We say that the Markov jump linear system (3) is pos¬ 
itive [18] if xo > 0 implies x(f) >0 for every t > 0 with 
probability one. For the system (3) to be positive, it is 
necessary and sufficient that all the matrices Ai, ..., A ,y are 
Metzler. For example, £ defined in (2) is a positive Markov 
jump linear system. The next stability condition is obtained 
in [18]. 

Proposition 2.4: Assume that the Markov jump linear 
system (3) is positive. Let If £ M. NxN be the infinitesimal 
generator of a. Then (3) is mean stable if and only if the 
matrix 

si = n T 0 /„ + diag(A 1; ... ,A n ) £ R nAfxnAf 

is Hurwitz stable. 

We also recall the next proposition, which gives a suffi¬ 
cient condition for the almost sure stability of not necessarily 
positive Markov jump linear systems. 

Proposition 2.5 ([19, Theorem 4.2]): Assume that a has 
the unique stationary distribution n. Then (3) is almost surely 
stable if E[ii(A k )] < 0. 

III. Disease Spread over Randomly Switched 
Networks with Independent Edges 

The aim of this section is to give an easy to apply sufficient 
condition for the stability of the switched A-intertwined 
model E under the assumption that the underlying randomly 
switched network has independent edges. 

A. Computational Difficulty 

Before presenting the main result we first observe the com¬ 
putational difficultly of the stability analysis. Theoretically, 
the next proposition completely solves the problem by giving 
a necessary and sufficient condition for mean stability via the 
eigenvalues of a matrix. 

Proposition 3.1: Let II £ E' Vx,v be the infinitesimal gen¬ 
erator of the Markov process a. Define 

sip = n T <g>/„ + jB diag(A^,... jAcfy). 

Then, E is mean stable if and only if 7] (.sip ) < S. 

Proof: The Markov jump linear system E is positive be¬ 
cause the matrix fiAg. — Si is Metzler for every i. Therefore, 
using Proposition 2.4, one can see that Tj(xip) < S implies 


that E is mean stable. Thus, by Lemma 2.2, we can show 
that E is also mean stable. 

On the other hand, assume that E is mean stable. Here 
we provide only a sketch of the proof. By Theorem 3.8 
in [20], we can show that E admits Lyapunov functions. 
Then it turns out that these Lyapunov functions also apply 
to E around a sufficiently small neighborhood of the origin. 
This is because, in such a neighborhood, the second order 
term [5PA<g a p of E can be ignored compared with its linear 
term. Therefore we can conclude that E is also mean stable. 
Hence, again by Proposition 2.4, we obtain rj (sip) < S. ■ 

Though Proposition 3.1 gives a practical characterization 
of stability when the size n of a switched network is 
effectively small, unfortunately, it cannot be easily used 
for real world networks with a large n. The size of the 
matrix .</3p equals nN, where A is the number of all the 
possible configurations of the switched network. This A can 
be very large when the network has many agents and can 
have various patterns of topology. In the extreme case when 
all the (") = n(n — l )/2 possible edges can be present or 
not independently, A equals the rapidly growing exponen¬ 
tial 2" 1 " l;/2 . By the same reason, Proposition 2.5 is often 
not easy to use in practice when n is large because it involves 
finding the matrix measure of a potentially huge number of 
matrices. We also remark that, in such a situation, it would 
be computationally too expensive to even estimate a stability 
condition by simulation. 

These difficulties motivate us to find easy to apply suffi¬ 
cient conditions for stability when the size of a network is 
large. 

B. Stability Condition for Switched Edge-independent Net¬ 
works 

One of the classical models of large-scale and random 
but static networks is Erdos-Renyi graphs [21], in which 
edges are assigned uniformly and independently for each 
pair of vertexes. Based on the fact that the distributions of 
the degrees of Erdos-Renyi graphs shows a large deviation 
from those of real world networks, recently Chung [11] 
proposed an improved version of the graphs by removing 
the uniformity constraint. 

Extending the above models, in this paper we study the 
disease spread over edge-independent random dynamical 
networks defined as follows. 

Definition 3.2: For distinct i.j £ { \.... .n} we let 
Aij = {Aij(t)} t > o denote the (i,y)-element of the matrix¬ 
valued stochastic process Ay a . We say that has inde¬ 
pendent edges if the processes {A (/ },>y are stochastically 
independent. 

Then our assumption on the randomly switched net¬ 
work !^ a can be stated as follows. 

Assumption 3.3: has independent edges and a has a 

unique stationary distribution. 

Remark 3.4: The existence of a unique stationary distri¬ 
bution is not very restrictive because a Markov process in 
general has a unique stationary distribution under a mild 
assumption of irreducibility and recurrence property. 


Let n denote the stationary distribution of a. Then 
is the random stationary graph with the random adjacency 
matrix Ay % . Define A = E\A<g n ] and 


A = 


max 

1 <i<n 


£Ay(l-Ay) 


;=i 


(4) 


Without loss of generality we assume A > 0, because other¬ 
wise A = 0 and therefore the graph process ^ equals the 
static network having the { 0 , l}-matrix A as its adjacency 
matrix. In this case E coincides with the AMntertwined 
model (1) for static networks and therefore E is stable if 
and only if A max (A) <8/j 3 as found in [5]. 

The next theorem gives an easy-to-use alternative of 
Proposition 3.1 and is the main result of this paper. 

Theorem 3.5: Define the function / on [0,°°) by 

/W= I + 2 n J exp(-^£ iI ). (5) 

Then E is almost surely stable if 

Amax (A) + min/(s) < 5/j8. ( 6 ) 

s>0 

Remark 3.6: We can understand the quantity mini >0 f(s) 
as the measure of uncertainty because the quantity increases 
with respect to A, which measures the variance of the 
stationary random graph ( S lt . 

We note that the size of the matrix A in the condition ( 6 ) 
is only n and is much smaller than that of sz7p in Propo¬ 
sition 3.1. Also A can be found as follows. Since has 
independent edges, each scalar-valued stochastic process A, ; 
is a {0, l}-valued time-homogeneous Markov process. Then 
its infinitesimal generator is of the form 

~Pij Pij 

qij -qij_ 


for some Pij,qtj > 0. The uniqueness of the stationary distri¬ 
bution of a yields that pij + qij > 0. Then we can easily show 
that the stationary distribution of Ajj, denoted by K t j, is given 
by ^({0})= /(pij + qij) and %({1 })= Pij/(Pij + qij). 
Therefore we obtain Ajj = p^/(pjj + qtj) when i^ j and 
An = 0 for every i. 

Moreover, the next proposition shows that the mini¬ 
mum min A >o/(,v) in ( 6 ) can be found by solving a convex 
program. 

Proposition 3.7: There exists 0 < so < 2A such that / is 
convex on [so, 00 ) and 

min/(s) = min(/(0),min/(s)) . (7) 


C. Proofs 

In this section we give the proofs of Theorem 3.5 and 
Proposition 3.7. For the proof of Theorem 3.5 we state two 
propositions. 

Proposition 3.8: If is [Amax (Ag^)] < 5/[3, then E is almost 
surely stable. 

Proof: Recall [13] that, if A is a real symmetric matrix, 
then ft (A) = Am ax (A). Thus we have ft(A^) = Am ax (A % ) 


with probability one. Now assume £ [A max (Ap^ )] < s/p. 
Then £[/i(/3A % - 5/)] = £[Vax(j3A % - Si)} < 0. There¬ 
fore, by Proposition 2.5, E is almost surely stable; i.e., 
||/>(f)|| —> 0 as t —> °° with probability one. Therefore, by 
using the equivalence of the norms ||-|| and ||-||i and also 
Lemma 2.2, we can show the almost sure stability of E. ■ 

To evaluate £[A max (A)] we will need the following result 
from the spectral theory of random graphs. 

Proposition 3.9 ([22, p. 7]): Let S? be a random undi¬ 
rected graph on the vertex set {1 ,..., n}, where two vertexes 
are adjacent in independently. Let A = £[A^] and define 
A by (4). Then, for every s > 0, 

P({A max (A) > A max (A) +s}) < 2n exp • ( 8 ) 

Now we prove Theorem 3.5. 

Proof of Theorem 3.5: Notice that the function / is 
well defined. Since has independent edges, any pair 
of vertexes are adjacent independently with others in the 
random graph ( S % . Therefore Proposition 3.9 applies to 
Let Q. denote the underlying probability space and define, 
for each s > 0, 

&S = { CO £ £2 . A max (A<g K ) i> Amax (A) T . 

If co € £2J, then A max (At^ ;l ) < A max (A) + s. On the other hand, 
if co € £2 S , then we have the trivial estimate A m ^ x (Ac^ n ) < n. 
Therefore, by (8) we obtain 

£[Amax(A % )] < P(f2j)(Amax(A) + s) +P{Cl s )n 

< Amax (A) + ^ + 2n 2 exp ^ ^ (9) 

= Amax (A ) +/(s). 

Since / is continuous and also f(s) diverges to +°° as s —> 
the minimum min J >o/(s) exists. Therefore, taking the min¬ 
imum with respect to s > 0 on the most right hand side of 
(9) proves £[A ma x(A^)] < Amax(A)+min s >o/(5). Hence, if 
(6) holds, then E is almost surely stable by Proposition 3.8. 


Then we give the proof of Proposition 3.7. 

Proof of Proposition 3.7: Let us only prove the convex¬ 
ity part, as the equation (7) follows from a straightforward 
argument. Throughout this proof we shall work with the 
translation of / defined by g(s ) = f(s — 3A) (s > 3A). We 
need to show that there exists 3A < sq < 5A such that g is 
convex on [so,°°). Define /r 1 ,/zt : [3A,°°) —> R by 


hi{s) 


Cl ( C 2 s 2 - c 3 + V2c 3 s) 
s 4 exp(c2« + C3S“ 1 ) 


ll2(s) = C2S 2 


C3-y/2^, 


where ci = 2 n 2 e 9A , C 2 = 3/2, and C 3 = 27A 2 /2. Then we 
can easily show that g" = h\li 2 - It is trivial to check that 
h\(s) > 0. Also, since /i 2 (3A) < 0 and h' 2 (s) > 0, there exists 
so > 3A such that /i 2 < 0 on [3A,so) and I 12 > 0 on (so, 00 ). 
Since /i(5A) = 24A 2 - Sv^A > (24 - 3v / lt)A 2 > 0, we can 
check so < 5A. The above argument yields g” < 0 on [3A, so), 
as desired. Then we can prove (7) by carefully investigating 
the derivatives g’ and g". The details are omitted. ■ 










IV. Examples 

In this section we apply Theorem 3.5 to the following 
important classes of dynamical graphs: randomly switched 
graphs with communities [23] and expected degrees [11]. 


A. Switched Graphs with Community Structure 

Let us consider a population grouped into the two com¬ 
munities Y\ = {l,...,ni} and Yi = {«i + l,...,«i +« 2 }- 
We assume that, for any pair of vertexes (i,j), the edge- 
process Ajj has the stationary distribution u, 7 that depends 
only on the communities to which i and j belong. Therefore, 
there exist 0 i, 02,0 £ [0,1] such that 


MOD 


9 e ijeVe, ^ = 1 , 2 , 

tj) otherwise. 


Then it follows that 


A = 


@\In[ 

0 


,«2 

^2^-«2,M2 ^42 


Therefore 

A = max((«i — 1 )0i (1 — 0i) +« 20(1 — 0), 

(n 2 - 1 ) 0 2 ( 1 - 62 ) + npj)(l -$)). 
Also we can show that 


■^max (-^) 


Hi 01 +H 202 + \/{n\ 6 \ -H 202 ) 2 + 4 «lH 20 2 
2 


— e 


for some e lying in between 0 ] and 02 - 

For example let n\ = 10 4 , «2 = 10 5 , 0i = 0.5, 02 = 0.3, 
and (j) = 0.1. Then one can compute A max (A) = 3.04- 10 4 . 
Also, by solving the convex program (7), we can easily find 
max 1 >o/(i) = 9.83- 10 2 , which is effectively smaller than 
Dnax (A). Then, by Theorem 3.5, E is almost surely stable if 
3.14 -10 4 < 8/15. We remark that, in this case, it is almost 
impossible to use Proposition 3.1 because the matrix has 
the dimension more than 10 10 ’. 


B. Switched Graphs with Expected Degrees 

In this section we study the special case when the expected 
adjacency matrix A has the following structure. 

Assumption 4.1: There exist ai,...,a n >0 such that 

Aij = ctjttj 

for every distinct pair (i.j). 

Epidemiologically, one can regard the constant a,- as 
the measure of the activity of the vertex i. Then we can 
understand Assumption 4.1 as stating that the frequency 
of the interaction between two agents i and j is solely 
determined by those activity measure. 

Also Assumption 4.1 can be supported by its connection 
to one of the well known models of random graphs. Let 
d £ M" be nonzero and nonnegative and let p = 1 /(£” =| d/). 
We say that an undirected random graph Y has expected 
degrees d [ 11 ] if edges are independently assigned to each 
pair of vertexes (i. j) with probability pdjdj. We can show 
the next proposition under Assumptions 3.3 and 4.1. The 
proof is straightforward and is hence omitted. 


Proposition 4.2: For each i let di = a, Y!j= 1 0 £/. The sta¬ 
tionary graph % has expected degrees d. 

The next theorem gives a sufficient condition for the 
almost sure stability of E in terms of the expected degrees 
of the stationary graph ( S K . 

Theorem 4.3: Let 

n n 

d = pYdf, Aj = max Y pdidj(p — djdj) 

and define the function / by (5) with A replaced by A,/. Then 
E is almost surely stable if 

d + min/(s) < 8//5. (10) 

^>0 

Moreover min i >o/(s) satisfies (7). 

Proof: It is straightforward to see that, for the random 
graph ( S K . the quantity A defined in (4) equals A f /. Moreover, 
since A = pdd J , we can show Amax^) = d. Thus the 
sufficient condition (10) immediately follows from (6). The 
proof of the latter claim is exactly the same as that of 
Proposition 3.7. ■ 

Example 4.4: One of the well-used model of the degree 
sequences is those with power law distributions [11], [22]. 
We say that a degree sequence d £ R" has the power law 
distribution with the power law exponent >2, maximum 
degree A > 0, and average degree d > 0 if 

di = c(i + io~ 1) - 1 /(P -1 ) 


for every i, where 

P ~ 2 Jn l W-V 


c = 


J 8 -1 


i 0 = n 


d(P~ 2 ) 
m(f5 — 1 ) 


/3-1 


For example let n = 10 7 , j5 = 2.2, A = 5 • 10 5 , and d = 10 3 . 
Then we have d = 3.15 • 10 4 . Also, by solving the convex 
program (7) we obtain min v ,o fif) = 1.97 • 10 3 . Therefore, by 
Theorem 4.3, E is almost surely stable if 3.35 • 10 4 < 5//J. 
We remark that the quantity min vx o/(.sj measuring uncer¬ 
tainty is relatively small compared with Amaxl/i). 

V. Disease Spread over Randomly Switched 
Networks with Weights 


In this section, we extend the result obtained in Section III 
to the networks modeled by weighted graphs. We start by 
giving necessary definitions. 

Definition 5.1: Let f be a finite set and let S be the 
set of all the distinct and unordered pairs of the elements 
of V. Also let w: $ —> [0,°°): (i.j) *-^w(i,j) be a function. 
We call the pair ('V. vv) a weighted undirected graph. An 
element of 'V (S’) is called a vertex (edge, respectively). For 
an edge e £ S, We call w(e) the weight of e. 

Let "V = {1,.The adjacency matrix App^ £ R' ,x ” 
of an weighted undirected graph (Y,w) is defined by 
[ A (y,w)]ij = w (hj) for i ^ j and l A (r,w)\u = 0 for every i. 
When no confusion arises we write A^-p^ as A w . We notice 
that A w is symmetric because (i,j) and (J, i) are the same 
unordered pairs. Then we call the differential equation 

-j- = (I5A W - 8l)p - f5PA w p 
at 








the A'-intertwined model of disease spread over ("V. w). 
Epidemiologic ally, the weight w(i,j) expresses the strength 
of the communication between the agents i and j, in the 
sense that instantaneous rate that disease transmits from i 
to j (provided i is infected and j is not) equals the prod¬ 
uct f3w(i,j). This in particular implies that, without loss of 
generality, we can normalize w as 

w(i,j) < 1 ( 11 ) 

by taking sufficiently large. 

Then, as in Section II-B, let us assume that the weight of 
the given network of agents changes over time according to 
a time-homogeneous Markov process, i.e., suppose that there 
exists a time-homogeneous Markov process <7 = {<7(f)};>o 
such that the weight of the graph at time t > 0 is given by the 
function w a ^ : £ —> [ 0 ,°°); {i.f) h. w a ^(i,j). Then we can 
model the disease spread over time-varying, weighted, and 
undirected graph & a := ("P. w a ) by the switched differential 
equation 

: —j~ = (/^w CT (,) — ^ ! )P — 

The almost sure stability of E u; is defined in the same 
way as Definition 2.1. Also, extending Definition 3.2, We 
say that ( £ a has independent edges if the stochastic pro¬ 
cesses j)} are independent. 

If a has a stationary distribution n, then we let Var(A H , ;I ) 
be the nxn real matrix obtained by taking the variances of 
the random matrix A Wn entry-wise. Since o is assumed to 
have finitely many states, each entry of A Wn is a distribution 
having finite support so that it is straightforward to find 
Var(A vlJI ). The next theorem extends Theorem 3.5 to the 
weighted and randomly switched networks. 

Theorem 5.2: Assume that ( £ a has independent edges and 
a has the unique stationary distribution 7 Z. Let A = E\A w f 
and also let A be the maximum row sum of Var(A WjI ). Define 
the function / by (5). If ( 6 ) holds, then L, v is almost surely 
stable. Moreover the minimum min s >o/(s) can be found by 
solving the convex program (7). 

Proof: From [22] we can easily check that Proposi¬ 
tion 3.9 holds for a weighted random graph as long as its 
weights are all less than one with probability one and we 
replace A by the maximum row sum of Var(A^). Therefore, 
under the normalization (11), we can apply Proposition 3.9 
to the stationary random graph ( S n and prove the theorem in 
the same way as Theorem 3.5. The details of the proof are 
omitted. ■ 

VI. Conclusion 

We studied the disease spread over randomly switched 
networks having stochastically independent edges. The dis¬ 
ease spread was modeled by a switched version of the 
V-intertwined model for static networks. Using the stability 
theory of Markov jump linear systems and the spectral theory 
of random matrices, we gave sufficient conditions for epi¬ 
demics dying out almost surely. We can check the proposed 
conditions by finding the maximum real eigenvalue of a 


matrix whose size equal the number of the agents and thus 
is highly efficient compared with another condition based 
solely on the stability theory. We also gave an extension to 
the case when networks are modeled by weighted graphs. 
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